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1.  INTRODUCTION 


A considerable  amount  of  effort  has  previously  been  devoted  to 
optimizing  the  performance  of  terrestrial  communication  systems  [1,2], 

The  objective  of  such  design  efforts  has  been  to  obtain  procedures  for 
optimizing  the  network  performance  with  regard  to  cost,  consumer 
satisfaction,  etc. 

With  the  growing  availability  of  commun Icatlon  satel 1 1 tes,  new 
optimization  procedures  for  the  hybrid  systems,  comprising  the  satellite 
and  terrestrial  communication  networks,  must  be  considered.  With  the 
growing  traffic  demands — and  the  poor  performance  of  some  of  the 
existing  terrestrial  system,  e.g.,  European  AUTOVON — new  algorithms  for 
Improving  network  performance  must  be  sought. 

The  hybrid  system,  like  the  terrestrial  systems,  could  be  used  in 
a fixed  network  configuration;  tnat  is,  the  satellite  capacity  available 
to  different  ground  stations  is  fixed  and  inflexible.  Such  an  operation 
would,  however,  ignore  the  most  important  feature  associated  with 
satellite  communications:  network  adaptability.  The  network  adaptability 

property  coupled  with  adaptive  routing  procedures  provides  extra 
flexibility  to  the  designer  to  improve  the  performance  of  the  hybrid 
system. 

One  aspect  of  commun icat ion  network  performance  optimization  deals 
with  the  routing  of  calls.  The  improvement  of  network  performance  for 
the  terrestrial  systems  through  routing  table  updating,  as  previously 
investigated  [3],  made  tacit  use  of  the  fact  that  the  network  config- 
uration is  fixed. 

With  the  satellite  system  the  network  topology  is  no  longer  fixed. 
Satellite  circuits  can  be  allocated  depending  on  the  demands  from  the 
ground  stations.  The  up-1  Ink  and  down-link  transmission  delay  prohibits 


any  alternate  (multiple  hop)  routing  on  the  satellite  system. 

Thus,  in  the  terrestrial  system  the  network  configuration  is 
fixed,  but  alternate  routing  is  possible.  In  the  satellite  system, 
the  network  topology  is  variable  but  only  direct  routing  is  allowed. 

In  the  hybrid  system,  therefore,  the  need  for  proper  use  of  these 
available  features  is  evident. 

This  report  documents  a preliminary  investigation  of  the  routing 
problem  in  hybrid  communication  networks.  A search  of  the  literature 
reveals  no  existing  work  on  such  networks. 

Section  2 deals  with  the  modelling  and  the  analysis  of  the  hybrid 
communication  networks. 

In  section  3,  the  optimization  problem  for  such  networks  is 
formu la  ted. 

Section  addresses  some  of  the  basic  questions  regarding  the 
effect  of  network  configuration  and  alternate  routing  on  trunk  group 
sizing.  A solution  to  the  problem  posed  in  section  3 is  also  presented 
in  this  section. 

Finally,  several  topics  of  future  research  are  given  in  section  5. 
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2 . MOPFLIMG,  ROUTING  STRATEGIES  AND  ANALYSIS  OF  HYBRID  COHflUN I CAT! ON  NETWORKS 
2. 1 Preliminary  Definitions 

A communication  network  can  be  modeled  by  a graph,  where  each  node 
in  the  graph  represents  a svj\  tchi ng  center  (exchange) , and  each  I ink 
represents  a trunk  group.  The  capaci ty  of  each  link  is  defined  in  terms 
of  the  number  of  trunks  in  the  link.  Each  t runk  represents  one  message 
channe 1 in  the  link. 

A network  with  n nodes  has  n x (n-1)  distinct  node  pairs.  The 
traffic  demands  for  all  node  pairs  can  be  described  by  a square  matrix 
of  orde:  n,  called  the  traffic  matrix.  The  diagonal  entries  of  this 
matrix  are  vacant  (have  no  physical  meaning).  Each  off-diagonal  element 
has  two  nodes  associated  with  it:  the  source  (S)  node  and  the  dest i na t ion 

(0)  node.  The  source  node  is  where  the  call  or i gi nate ; and  the  destination 
node  is  where  the  calls  terminate  (or  are  destined  to).  The  elements  of 
the  traffic  matrix  may  be  in  units  of  Erlangs  or  hundred-call  seconds 
(abbreviated  CCS),  one  CCS  being  equal  to  l/36th  of  an  Erlang, 

The  class  of  communication  networks  dealt  with  in  this  report 
belong  to  the  circuit-switched  type.  In  circuit-switched  networks,  a 
path  is  set  up  between  the  source  and  the  destination  pair  prior  to  the 
starting  of  message  transmission.' 

The  nodes  are  assumed  to  be  perfectly  reliable  (blocking  probability 
= O).  The  links,  however,  are  nwdeled  on  the  basis  of  their  capacity, 
and  their  availability  is,  therefore,  assumed  to  be  probabilistic.  Based 
on  this  model,  the  link  grade  of  service  is  defined  as  the  calls  blocked 
by  a link  divided  by  the  calls  attempted  on  the  link. 

'por  example,  in  the  SPADE  system,  the  path  (channel)  is  set  up 
prior  to  message  transmission. 
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Due  to  the  probablMstic  nature  of  the  avaMability  of  links,  the  ! 

t 

actual  calls  completed  between  a source-to-dest Inat ion  pair  will  be  1 

less  than  the  originating  calls.  To  this  end,  we  define  the  node-to- 

node  grade  of  service  (abbreviated  NNGOS)  as  the  lost  load  between  a 

source-+dest inat ion  pair  divided  by  the  total  load  destined  between  the 

sou rce'^est inat ion  pair.  Two  other  quantities  that  are  often  used  in 

describing  the  network  performance  are:  (1)  the  network  grade  of  service 

is  defined  as  the  total  traffic  lost  in  the  network  divided  by  the  total 

2 

traffic  originating  in  the  network,  and  (2)  the  node  grade  of  service 
is  defined  as  the  total  traffic  lost  at  a node  divided  by  the  total 
traffic  originating  at  the  node. 

2.2  Routing  Strategies 

The  routing  strategy  employed  in  a network  can  be  completely  des- 
cr i bed  by  a routing  table  and  a cal  1 control  rule. 

The  routing  table  for  an  n-node  network  can  be  expressed  by  a 
matrix  of  order  n.  This  matrix  is  used  to  dictate  the  "preference 
scheme"  for  routing  calls.  The  order  of  the  entries  in  the  routing 
table  indicates  the  preference  scheme. 

The  need  to  define  a call  control  rule  arises  due  to  ambiguities 
that  may  arise  in  the  routing  table  spec  if icat ion  [3].  A number  of  call 
control  rules  exist  in  various  communication  networks  [A].  In  this 
report,  we  shall  deal  with  the  originating  office  control  rule.  The 
originating  office  control  (abbreviated  OOC)  rule  may  be  stated  as  follows: 

2 

Note  that  this  quantity  is  due  to  the  probabilistic  nature  of  the 
availability  of  the  1 i nks  and  not  the  nodes  (as  assumed  earlier). 
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"If  node  I is  the  originating  node  for  a call  destined  to  node  D, 
then  the  call  is  routed  to  some  adjacent  node  specified  by  the  from- 1 - 
to-D  block  in  the  routing  table.  If  I is  a tandem  node,  then  a call 
reaching  node  I is  routed  only  to  the  node  appearing  as  the  f i rst 
choice  in  the  from-l-to-D  block  in  the  routing  table.  If  the  link 
associated  with  this  choice  is  unavailable,  then  the  originating 
office  (the  source  node)  is  informed  of  this  condition  and  the  call 
is  attempted  on  the  next  link  going  from  the  source  node.  If  no 
such  link  exists  at  the  originating  office,  the  call  is  lost." 

In  some  networks,  the  number  of  links  incident  to  a node  may 
be  small.  In  order  to  give  added  flexibility  to  the  originating  office 
to  route  calls,  the  OOC  rule  is  modified  to  allow  for  spill  forward 
action.  According  to  the  originating  office  control  with  spill  forward 
rule,  the  originating  office  may  pass  (spill)  the  control  of  a call  to 
a tandem  switch,  which  then  assumes  the  role  of  the  originating  office. 

In  general,  the  number  of  paths  available  to  route  calls  increases  with 
the  spill  forward  action. 

The  process  of  routing  can  be  depicted  with  the  aid  of  a route  tree, 
which  essentially  shows  all  the  paths  defined  by  the  routing  table  and 
the  call  control  rule.  In  order  to  completely  define  the  routing 
procedure,  the  concept  of  augmented  route  tree  is  used  [3].  The  aug- 
mented route  tree  can  be  obtained  from  a route  tree  in  the  following 
manner: 

Add  one  directed  branch  below  ail  the  branches  leaving  an 
originating  node  or  the  spill  node  in  the  route  tree,  and  label  this 


I 

) 


j 
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terminal  node  as  L. 
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The  route  tree  of  Fig.  2.1(a)  depicts  the  routing  of  calls  from 
node  to  node  B.  The  augmented  route  tree  of  Fig.  2.1(b)  depicts  the 
routing  as  well  as  the  loss  of  calls  for  the  same  node  pair. 


B B 


(a)  Route  Tree 


(b)  Augmented  Route  Tree 


Figure  2.1 


2,3  Hybrid  Communication  System  Network  Structure 

The  definitions  Introduced  in  section  2.1  are  sufficient  to  model 
a terrestrial  communication  systems.  The  same  modelling  scheme  can  be 
used  to  represent  a hybrid  communication  system.  However,  if  a demand 
assignment  scheme  is  to  be  implemented  on  a digital  computer,  the  need 
to  distinguish  the  satellite  portion  of  the  hybrid  system  is  evident. 

In  this  section,  we  will  present  a network  model  that  can  be  used  to 
represent  any  hybrid  communication  system. 

In  the  terrestrial  communication  system,  the  nodes  represent  a 
switching  center  of  the  network,  and  a link  is  used  to  connect  two 
switching  centers.  In  such  networks,  the  link  capacity  and  the  location 
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of  the  switching  centers  are  fixed.  Thus,  the  network  has  a fixed 
conf igurat ion. 

In  the  satellite  communication  system,  the  nodes  represent  a 
transit  center  (TC)  or  a ground  station  (GS).  Due  to  the  large 
expense  involved  in  obtaining  and  maintaining  a ground  station,  not 
all  the  countries  have  such  facilities.  Such  countries  can  access 
the  satellite  by  routing  the  call  to  the  nearest  ground  station.  Thus 
the  links  in  such  networks  represent  connections  between  two  transit 
centers,  two  ground  stations,  or  a transit  center  and  a ground  station. 
The  capacities  of  the  links  connecting  two  transit  centers  or  transit 
center  and  ground  station  are,  in  general,  fixed.  The  capacity  of  the 
link  whose  both  ends  terminate  at  a ground  station  are  variable. 

The  hybrid  communication  system,  which  is  a combination  of  the 
terrestrial  and  the  satellite  system,  can  be  represented  by  the  block 
diagram  of  Fig.  2.2. 


r- 
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I 
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1 
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TERRESTRIAL  SATELLITE 
LINKING  NETWORK 
(fixed  configuration) 

Figure  2.2.  A Hybrid  Communication  System 
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■}  In  modeling  the  network  under  study,  the  following  links  may  be 

1 

I considered:^ 

1 

I 1.  Pre-assigned  Trunks:  Such  trunks  may  connect  a transit 

^ center  (or  switching  center)  to  another  transit  center  or 

J 

' a ground  station.  These  links  operate  between  weli-defined 

i 

j fixed  points,  and  their  capacities  are  fixed. 

2.  Fully  Variable  Trunks:  Both  the  terminal  nodes  of  such 

, trunks  are  the  nodes  representing  ground  stations.  Their 

i 

I capacities  as  well  as  their  terminal  nodes  are  variable. 

j 

j The  total  number  of  fully  variable  trunks  in  the  network 

■ are  fixed. 

With  the  two  types  of  trunks  and  two  types  of  nodes  (TC  and  GS) , 
a hybrid  communication  network  can  be  represented  using  four  different 
symbols.  An  example  of  such  a network  is  shown  in  Fig.  2.3.  The 
circular  nodes  represent  the  transit  centers  (or  switching  centers) 
and  the  triangular  nodes  represent  the  ground  stations.  The  solid 
lines  represent  the  pre-assigned  trunks;  while  the  broken  (dashed) 
lines  represent  the  fully  variable  trunks.  A location  having  both 
the  transit  center  and  the  ground  station  in  close  proximity  is 
represented  by  a circle  inserted  inside  the  triangle. 

The  routing  table  for  the  hybrid  network  is  constructed  in  a 
manner  similar  to  that  for  terrestrial  networks.  A country  without 
a ground  station  facility  can  spill  its  calls  to  the  closest  earth 
station.  For  example,  in  Fig.  2.3,  the  switching  center  represented 

^If  the  economic  aspects  are  also  considered,  then  the  need 
for  variable  destination  links  must  also  be  considered.  Such  links 
arise  when  the  1 inks  are  leased  by  a country.  However,  our  study 
will  be  concerned  with  the  two  types  mentioned  above. 


Figure  2.3.  A Hybrid  Communication  Network 
by  node  number  1 can  spill  all  of  its  calls  to  the  ground  station 

numbered  b.  By  spilling  the  call,  node  1 can  also  share  the  "benefits"  ' 

of  alternate  routing. 

2.4  Analysis  Method 

Most  of  the  optimizing  criteria  for  communication  networks  are 
related  to  improvement  in  the  network  performance,  reduction  in  cost, 

etc.  In  this  study,  the  optimization  criteria  is  the  improvement  of  i 

1 

the  network  performance.  An  in  depth  discussion  of  the  problem  is  ^ 

presented  in  the  next  section. 

In  order  to  improve  network  performance,  a method  of  some  “how 
estimating  the  performance  is  needed.  Two  kinds  of  methods  are  generally 
used  for  such  studies:  simulation  and  analytic  [5].  In  the  former 

method,  artificial  traffic  is  generated  on  the  computer,  and  an  estimate 
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of  the  system  behavior  is  obtained.  Simulation  methods  may  look 
attractive  because  of  the  inherent  flexibility.  However,  when  large 
scale  networks  are  simulated  on  the  computer,  the  advantage  due  to 
the  flexibility  may  be  abated  by  the  high  cost  involved  in  the 
sijnulat  ion. 

In  the  analytical  methods,  the  only  disadvantage  is  due  to  the 
approximations  made  in  deriving  the  method.  Most  of  the  studies  have 
shown  that  the  results  obtained  from  analytical  methods  are  in  good 
agreement  with  those  of  the  simulation  methods.  The  advantage  of 
analytical  methods  is  the  small  computation  tinr.e  required  for 
estimating  the  network  performance. 

In  the  demand  assignment  mode  of  operation,  it  is  important 
that  a fast  estimate  of  the  network  performance  be  obtained,  so  that 
the  fully  variable  trunks  can  be  re-allocated  to  meet  the  changing 
traffic  demands.  Consequently,  for  the  problem  to  be  presented  in 
the  next  section,  analytical  techniques  are  preferred. 

The  analytical  method  used  for  obtaining  the  network  performance 
has  been  adapted  from  [3].  An  overview  of  this  so-called  single- 
moment method  will  be  presented  here. 

In  deriving  the  singlf*-moment  method,  the  following  assumptions 
are  made: 

1.  Call  Arrival  is  a Poisson  process 

2.  Call  holding  time  has  a negative  exponential  distribution 

3.  Link  blocking  probabilites  are  statistically  independent 

U.  Nodes  are  non-blocking  (has  been  stated  earlier) 

5.  Blocked  calls  are  cleared  and  do  not  return 

6.  The  network  is  in  statistical  equilibrium 
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i 
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7.  Call  set-up  time  is  negligible 

8.  Call  arrival  on  any  1 i nk  is  a Poisson  process 

By  making  the  assumption  that  the  overflow  traffic  is  also  Poisson 
distributed  (assumption  8),  we  obviate  the  need  to  use  the  dual-moment 
method  [5].  By  virture  of  this  assumption,  the  mean  and  the  variance 
are  equal,  and  the  process  is  characterized  only  by  the  mean — hence, 
the  name  "single-moment  method," 

The  following  steps  are  used  to  determine  all  the  node-to-node 
grade  of  services  by  the  single  moment  method: 

Step  1 : From  the  given  network  configuration,  the  call  control 

rule  (OOC  with  spill,  in  our  case),  and  the  routing  table,  find  the 
augmented  route  tree  for  every  node  pair. 

Step  2:  From  the  traffic  matrix  and  the  augmented  route  trees, 

find  the  link  blocking  probabilities  using  the  Erlang's  loss  formula  [6]: 

N. 

a.'  /N.! 

I I 

y = ^ i = 1,  2,  3,  ...  ^ (2.1) 

'^i  k 

I (a^/k!) 

k=0 


where  a.  and  Nj  are  the  offered  loads  (in  Erlangs)  and  the  number  of 
trunks  in  the  ith  link,  respectively,  and  I is  the  total  number  of 
links  in  the  networks. 

Step  3:  From  the  link  blocking  probabilities  found  in  Step  2, 


i 

t 
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find  the  probability  of  each  path  in  the  augmented  route  tree  being 
used  to  complete  a call. 


From  the  above  steps  all  information  pertaining  to  the  network 
performance  can  be  obtained.  This  information  may  inciude:  node  grade 

of  service,  network  grade  of  service,  trunk  group  offered  loads,  load 
carried  by  each  path  between  every  node  pair,  etc.  For  a detailed 

i 

! description  of  this  algorithm,  refer  to  [3]. 

( 


I 


This  step  will  be  discussed  in  more  detail  in  section  *4. 
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3.  OPTmiZATION  PROBLEM 

The  design  of  circuit-switched  networks  involves  the  following 
basic  steps: 

(a)  Determination  of  network  configuration  (or  topology)  to  provide 
communication  services  to  the  users. 

(fa)  Determination  of  trunk  group  (iink)  sizes  and  routing  table  to 
provide  satisfactory  network  performance. 

In  the  design  of  terrestrial  communication  networks,  previous  approaches 
have  been  based  on  the  assumption  that  the  network  configuration  is  fixed. 

Once  the  network  topology  and  trunk  group  sizes  have  been  decided,  there  is 
little  freedom  to  improve  the  network  performance  by  altering  the  routing 
table.  Adaptive  routing  techniques  are  available  to  improve  the  network 
performance  [3] . But  when  the  traffic  load  becomes  too  high,  these  techniques 
are  inadequate  to  meet  the  demands.  The  optimization  criterion  used  in 
terrestrial  networks  can  be  reduction  in  cost,  improvement  of  the  network 
performance,  or  the  improvement  in  the  point-to-point  performance. 

A sateliite  communication  system  can  be  used  in  a changing  network 
configuration  mode.  In  fact,  to  take  advantage  of  the  inherent  flexibility 
of  such  systems,  network  adaptation  should  be  an  essential  control  feature. 
Along  with  the  network  adaptability  feature  is  the  freedom  to  vary  the  size 
of  the  variable  trunks.  The  only  restriction  is  that  the  total  sum  of  these 
variable  trunks  must  not  exceed  the  satellite  capacity. 

In  the  hybrid  communication  network,  therefore,  there  are  several  options 
available  to  the  designer  to  meet  a specified  network  performance.  Adaptive 
routing  techniques,  network  adaptability  and  the  allowed  changes  in  trunk 
group  sizes  can  be  used  to  meet  an  optimization  criterion. 
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The  optimization  criterion  that  we  have  chosen  in  this  report  is  totrltirlrc 
node-to-node  grade  of  service.  The  otjective  is  to  design  a network  for  a 
given  traffic  matrix  such  that  ali  the  node-to-node  grade  of  services  are 
beiow  a prescribed  vaiue  (abbreviated  BPMAX--maximum  aliowed  blocking 
probabi 1 i ty ) . 

To  give  a formal  statement  of  the  problem,  we  make  the  following  definitions; 

G = number  of  terrestrial  communication  nodes 

S = number  of  satellite  ground  stations 

T = [t  .]  = traffic  demand  matrix,  where  the  element  t..  specifies  the 

- T I j T U 

demand  from  node  i to  node  j at  time  T. 

C = [c.]  = link  capacity  vector,  where  the  C.  element  denotes  the  number 

^ T I T I 

of  trunks  in  the  ith  link  at  time  T. 

BPMAX  = maximum  allowable  blocking  probability  between  each  node  pair. 

= routing  table  at  timer. 

B = [b..]  = node-to-node  grade  of  service  matrix,  where  b. . element 
~ ij  ij 

denotes  the  NNGOS  of  i-to-j  node  pair. 

We  further  assume  that  some  of  the  elements  of  C are  fixed.  These 

~T 

correspond  to  the  links  between  a switching  c'^r.i.er  and  ground  station  or 
between  two  switching  centers. 

Our  objective  is  to  find  a routing  table  and  a link  capacity  vector  such 
that  all  the  node-to-node  grade  of  services  are  below  a prescribed  BPMAX. 
Furthermore,  in  sat i sf y ing th i s criterion,  the  number  of  fully  variable  trunks 
used  should  be  minimized.  Mathematically, 


C = {[c.]  lb.,  ^ BPMAX,  min  [Ec .] } 

T J I K J 

R = { [r.,  ] |b,,  < BPMAX} 

T ik  ' ik 


j = 1,  2,  3,  ...,  H 
i = 1,  2,  3»  ...,  G+S 
k=  1,  2,  3,  ...,  G+S 


Other  optimization  criteria  do  exist  [7]  but  we  will  be  concerned  with 


the  one  mentioned  above 


iround  stations.  The  up-link  and  down-link  transmission  delay  prohibits 


To  see  the  complexity  of  this  problem,  let  us  quantify  some  of  the 

available  options:  For  the  S node  satellite  system  there  can  be  a total  of 

S { S“  1 ' 

— links.  The  number  of  possible  link  combinations  is  identical  to  the 

number  of  possible  network  configurations.  This  upper  bound  can  be  easily 
derived  as: 


(?)  = 2^-1 


- S(S-l) 
where  N = — 


Each  of  the  network  configurations  can  be  further  modified  by  the  large 
number  of  possible  trunk  assignments.  Obviously,  there  is  an  astronomical 
number  of  the  combinations  of  conceivable  network  topologies  and  trunk 
group  sizing. 

The  traffic  matrix  and  the  prescribed  BPMAX  does  reduce  the  number  of 
the  possible  combinations.  Other  constraints  such  as  satellite/ground 
station  visibility,  power  sharing,  uplink  and  downlink  losses,  etc.  further 
reduce  the  upper  bound  on  all  the  possible  combinations.  Despite  all  of  these 
reductions  the  number  of  available  options  is  still  astronomical. 

Sometimes,  even  if  ^ solution  to  the  aforementioned  problem  is  found, 
there  are  other  problems  that  may  arise.  For  example,  it  is  possible  for  a 
calculated  route  to  use  the  same  satellite  more  than  once.  Such  a route  is, 
of  course,  undesirable  because  of  the  double  transmission  time  delay  involved. 
Should  such  a problem  arise,  the  need  to  re-al locate  link  capacities  and  obtain 
a new  routing  table  is  evident. 

Finally,  the  algorithm  for  solving  this  network  design  problem  must  be 
fast.  If  the  demand  requests  are  monitored  every  seconds,  the  algorithm 
should  be  able  to  make  the  routing  table  and  trunk  assignments  in  less  than 
\t  seconds  so  that  it  can  track  the  demand  assignment  monitor.  Such  an 
algorithm  could  then  be  used  in  real-time  computing. 
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li.  ALGORITHM  FOR  TRUNK  ASSIGNMENT  AND  ROUTING  TABLE  GENERATION 
k, 1 Introduction 

In  this  section  we  will  present  a heuristic  algorithm  for  obtaining  the 
branch  capacity  vector  and  the  routing  table,  for  a given  traffic  demand, 
that  will  meet  the  prescribed  BPMAX  specification.  The  algorithm  also 
attempts  to  minimize  the  number  of  fully  variable  trunks  used  in  the  network. 

The  next  two  sub-sections  are  devoted  to  showing  the  effect  of  network 
topology  and  alternate  routing  on  trunk  group  sizing.  The  ideas  gathered 
in  these  sub-sections  will  aid  in  the  development  of  the  algorithm  to  be 
presented  in  section 

i*.2  Effect  of  Network  Topology  on  the  Link  Capacity  Vector 

In  terrestrial  communication  networks  it  is  economically  unfeasible 
to  connect  a trunk  group  between  every  pair  of  switching  centers  in  the 
network.  Some  of  the  calls  have  to  be  either  spilled  to  another  office, 
or  to  be  routed  through  a tandem  office. 

N 

In  the  satellite  system,  one  of  the  (2  -1)  possible  network  con- 

.C 

figurations  is  a complete  gr^ph;  that  is,  there  Is  a link  between  every 
pair  of  ground  stations.  Consequently,  a portion  of  the  hybrid  communication 
network  that  we  are  concenred  with  can  be  represented  by  a complete  graph. 

An  obvious  question,  which  must  be  answered  before  any  decision  on 
the  trunk  group  sizing  algorithm  is  made,  is  the  following:  Is  it 

advantageous  to  operate  a network  in  a complete  graph — as  opposed  to  an 
Incomplete  graph — configuration  when  such  a choice  exists? 

See  Eqn.  (3.2) . 


A pure  graph-theoretic  point  of  view  reveais  that,  for  a n node 
network,  the  complete  graph  wi i 1 be  strongiy  connected  as  compared  to  an 
incompiete  graph.  Accordingiy,  it  will  be  less  vulnerable  or  more  damage 
resistant  [8]. 

At  this  time,  it  is  not  possible  to  give  a mathematical  proof,  based 
on  the  single-moment  method,'  to  show  that  the  total  number  of  trunks 
required  to  achieve  a prescribed  BPMAX  is,  in  general,  smaller  for  complete 
graphs  than  incomplete  graphs.  All  of  the  network  examples  analysed  using 
STARTUP  [3]  have  validated  this  conjecture.  Consider,  for  example,  the 
five-node  network  of  Fig.  4.1  (a).  Using  the  symbols  Introduced  in  section 
2.3,  the  pre-assigned  links  are  represented  by  solid  lines,  and  the  fully 
variable  links  by  broken  lines.  The  satellite  ground  stations  and  switching 
centers  of  the  terrestrial  network  are  grouped  together.  Notice  that  there 
are  no  terrestrial  links  between  node  pairs  (1,3)  and  (2,5). 

The  BPMAX  specification  is  0.025 — that  is,  all  the  node-to-node 
grade  of  services  must  be  less  than  2.5%.  The  traffic  matrix  for  this 
network  is  given  in  Table  4.1.  Assuming  that  the  combined  capacity 
of  the  terrestrial  network  is  not  enoggh  to  handle  the  traffic  demand, 
we  must  assign  satellite  trunks  (circuits)  to  fulfill  the  prescribed 
BPMAX  requirement. 

Two  of  the  several  network  configurations  possible  are  depicted 
in  Fig.  4.1(a)  and  (b).  In  Fig.  4.1(a)  we  have  assigned  the  fully 
variable  trunks  such  that  a complete  graph  is  formed.  In  Fig.  4.1(b) 
we  assign  trunks  to  increase  the  capacity  of  the  already  existing 
terrestrial  trunk  groups. 


This  is  a nonlinear  approach  to  the  performance  analysis;  based  on 
the  flow  analysis  approach,  one  could  show  this  conjecture  using  the  min 
cut-max  flow  argument. 
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Clearly,  there  are  several  different  combined  branch  capacity 
vectors  that  can  be  chosen  to  meet  the  BPMAX  specification.  One  of 
such  vectors  is  given  in  Table  l4.2(b).  Based  on  the  traffic  matrix  and 
the  trunk  group  sizes,  weobtain  the  routing  table  depicted  In  Table  4.2(a) 
such  that  all  the  NNGOS  are  less  than  0.025.  The  network  performance 
using  STARTUP  is  shown  in  Fig.  4.2. 
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Table  4.2  Routing  Table  and  Trunk  Group  Information  for  Fig.  4.1(a) 


(i)  Routing  Table 


Link  Number 

Terminal  Nodes 

Number  of  Trunks 

1 

1-2 

31 

2 

1-4 

44 

3 

1-5 

36 

4 

1-3 

41 

5 

2-3 

39 

6 

2-4 

45 

7 

2-5 

33 

8 

3-4 

28 

9 

3-5 

48 

10 

4-5 

38 

(b)  Trunk  Group  Information 

2 

By  combined  we  mean  the  sum  of  the  pre-assigned  and  fully  variable 
trunks  present  in  a link. 


L 


PCTUhL  TkUHK  GROUP  PERFORMArCE  (BASED  DM  ORIGINAL  ROUTING  TABLE) 


NETUORK  GRhDE  of  service  IS  .‘1177'=; 

UEIGHTED  rCTUOPF  GRADE  DP  SERVICE  IS  .01775 
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From  Figs.  4,2  and  4.3,  we  note  that  the  BPMAX  criterion  Is  met  in 
both  the  cases.  The  total  number  of  trunks  needed  for  the  network  config- 
urations of  Figs.  4.1(a)  and  (b)  is  3^3  and  403,  respectively.  Thus, 
there  is  a saving  of  trunks  if  we  use  the  complete  graph  configuration; 
furthermore,  the  network  grade  of  service  is  smaller  for  this  case. 

With  this  small  network  example,  we  have  illustrated  the  advantage 
accrued  in  choosing  a complete  graph  topology  over  another  topology.  As 
mentioned  earliei — even  for  this  simple  network — there  are  several  choices 

, for  the  branch  capaci  ty  vectors  for  a given  configuration.  In  all  of  the 

I 

! examples  investigated,  our  results  Indicate  that  the  complete  graph 

configuration  always  uses  fewer  trunks  than  another  configuration  to  meet 
the  same  BPMAX  specification.  A derivation  of  this  conjecture  should  be 
a part  of  any  future  research. 

4. 3 Effect  of  Alternate  Routing  on  the  Link  Capacity  Vector  ' 

In  this  section,  we  will  explore  the  effect  of  alternate  routing 
on  trunk  group  sizing.  If  we  were  concenred  only  with  the  satellite 
system,  such  an  investigation  would  be  unwarranted;  for,  the  Inherent 
delay  in  the  uplink  and  downlink  transmission  makes  the  use  of 

alternate  routing  strategy  undesirable.  However,  in  the  hybrid  system  \ 

alternate  routing  is  feasible.  We  can  assume  that  a call  is  routed  in  ^ 

a manner  such  that  is  uses  the  satel  1 i te  "trunk"  at  most  once  in  its 

propagation  from  the  source  to  the  destination.  Consequently,  it  is 

meaningful  to  investigate  the  behavior  of  NNGOS  for  different  trunk 

group  sizes  and  routing  strategies. 


2 <4 

The  relationship  between  the  offered  traffic  and  node-to-node  grade 
of  services  for  fixed  trunk  group  sizes  and  different  routing  schemes  has 
been  reported  in  [if],  Grandjean's  results  show  that  the  "nondirect"  routing 
with  more  alternate  paths  is  more  sensitive  to  overload  than  alternate 
routing  with  fewer  paths.  Direct  routing  has  the  least  sensitivity  with 
regard  to  overload. 

In  our  problem,  we  are  not  only  concerned  with  selecting  a routing 
strategy,  but  also  minimizing  the  number  of  trunks  required  to  meet  a 
prescribed  BPMAX  value.  An  obvious  question  that  arises  is:  For  a pre- 

scribed traffic  matrix,  what  is  the  relationship  between  the  node-to-node 
grade  of  services  and  varying  link  sizes  for  different  routing  schemes? 

We  will  show  this  reiationship  using  an  idealized  network  example. 
Consider  the  five-node  network  of  Fig.  4.^.  The  network  is  completely 
symmetrical  with  regard  to  traffic  and  routing  strategies;  that  is,  the 
traffic  between  each  node  pair  is  equal  and  every  link  has  the  same 
offered  and  carried  loads.  The  number  of  trunks  assigned  to  each  link 
is  the  same. 

2 


Figure  4. A.  A five-node  symmetrical  network 
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35 


3 


39 


Fig.  '4.5  NNGOS  vs. 

(Traffic 


NUMBER  OF  TRUNKS 

Nunber  of  Trunks  for  the  network  of  Fig.  k.ii 
= 15  Erlangs) 
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We  investigate  the  relationship  for  direct  routing  and  alternate 
routing  with  one  jnd  two  alternate  f>aths.  Because  of  the  symmetry  in  the 
traffic,  topology,  and  routing,  all  the  link  blocking  probabilities  and  the 
node-to-node  grades  of  service  for  every  node  pair  are  the  same. 

Fig.  h.S  depicts  the  relationship  between  the  node-to-node  grades  of 
service  and  varying  trunk  group  sizes  for  zero,  one,  and  two  alternate 
paths  between  every  node  pair.  It  is  interesting  to  note  some  of  the 
salient  features  of  Fig.  4.5: 

(a)  It  is  not  always  preferrable  to  use  alternate  routing  in  complete 
graphs.  The  routing  strategy  should  be  chosen  from  the  prescribed 
BPMAX  value  and  the  total  number  of  fully  variable  trunks  available. 
At  times  it  may  so  happen  that  the  traffic  demand  may  become 
excessively  high  so  that  even  with  all  the  fully  variable  trunks 
allocated  the  BPMAX  specification  is  not  satisfied.  The  network 
performance  can  still  be  improved  by  a proper  choice  of  the 
routing  strategy. 

(b)  The  two  alternate  route  strategy  is  more  sensitive  to  decreasing 
trunk  group  sizes  than  one  alternate  route  or  direct  strategies. 

The  higher  sensitivity  of  alternate  routing  due  to  decreasing 
link  capacities  is  generally  attributed  to  link  "congestion.” 

Even  though  we  have  used  an  idealized  example  to  illustrate  the 
dependence  of  the  node-to-node  grades  of  service  on  the  trunk  group  sizes 
for  different  routing  strategies,  the  results  obtained  are  applicable  to 
the  "real  world"  hybrid  communication  networks,  too.  As  will  become 
apparent  in  the  next  section,  our  algorithm  for  demand  assignment  attempts 
to  distribute  the  traffic  in  a manner  so  that  the  variation  in  the  offered 


loads  to  the  links  is  minimal. 
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A. 3 An  Algorithm 

A search  of  the  literature  reveals  no  existing  worl;  dealing  with 
the  solution  to  the  problem  formulated  in  section  3""that  of  determining 
the  trunk  group  sizes  and  routing  table  generation  s imul taneous 1 y.  Previous 
works  [3,9]  have  treated  these  as  disjoint  problems.  Reference  [3]  gives 
an  algorithm  for  obtaining  a routing  table  that  satisfies  a prescribed 
BPMAX  specification,  assuming  that  the  network  configuration  and  trunk  group 
sizes  are  fixed.  On  the  other  hand,  [9]  uses  an  iterative  scheme  to  obtain 
the  branch  capacity  vector  for  a pre-spec i f ied  routing  table. 

In  this  section,  we  will  present  an  algorithm  that  will  generate  a 
branch  capacity  vector  and  routing  table  for  a network  with  adaptable 
conf igurat ion,  such  that  all  the  NNGOS  are  below  a prescribed  BPMAX  value. 

Of  course,  if  the  traffic  demands  get  excessively  high  in  relation  to  the 
network  cal  1 -carry i ng  capacity,  this  algorithm  fails  to  meet  the  BPMAX 
cr i ter  ion. 

The  researcli  reported  herein  is  a preliminary  investigation  of  the 
routing  problem  in  hybrid  communication  systems.  Consequently,  no 
rigorous  derivations  of  the  concepts  used  are  available  at  this  point. 

We  shall,  however,  present  heuristic  arguments  to  justify  the  steps 
involved. 

Our  aim  is  to  determine  an  algorithmic  way  of  designing  networks  with 
"semi-adaptable"  configurations  such  that  ail  the  NNGOS  are  below  an  upper 
bound.  In  order  to  facilitate  the  understanding  of  the  various  steps,  and 
their  order,  in  the  algorithm,  we  will  show  how  the  NNGOS  are  calculated  and 
their  dependence  on  the  link  grade  of  services,  traffic  matrix,  and  the 


routing  table 
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A. 3.1  Factors  that  Influence  NNGOS  Values 

The  complexity  Inherent  in  the  hybrid  communication  network  routing 
problem  was  presented  in  section  3.  This  section  reveals  some  additional 
difficulties  involved  in  the  solution  to  the  problem  we  have  formulated. 

A systematic  approach  can  be  devised  after  the  basic  complexities  in  the 
problem  are  understood. 

The  network  performance  (determination  of  the  NNGOS)  based  on  the 
single-moment  method  can  be  obtained  provided  we  have  the  following 
i nformat ion: 

1.  Network  Topology 

(a)  link-node  incidence  description  (b)  link  sizes 

2.  Routing  Table 

3.  Call  Control  Rule 
A.  Traffic  Matrix 

In  our  case,  the  desired  network  performance,  (3)  and  (A)  are 

known;  the  problem  is  to  find  (1)  and  (2).  Clearly,  we  note  that  we 

have  an  extra  unknown  in  our  problem.  Thus,  an  obvious  question  is: 

knowing  the  desired  network  performance  (specified  by  the  BPMAX  requirement), 

(3)  and  (A)  above,  is  It  possible  to  find  (1)  and  (2)?  This  question  can 

be  answered  after  we  develop  the  method  for  determining  the  NNGOS  for  all 

H • ^ 

node  pairs:  j 

We  first  define  the  following  notation: 

vector  of  link  blocking  probabilities  ^ 

= vector  ofllnkreliabilities  j 

i 

a^'  = vector  of  link  carried  loads  j 

vector  of  link  offered  loads  | 

A . . . ! 

c = branch  capacity  vector  , 


I 
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In  order  to  determine  the  NNGOS  for  all  node  pairs,  the  ^vector, 

the  call  control  rule,  and  the  routing  table  must  be  known.  The  ^vector 

is  found  by  using  the  Erlang  loss  formula  [6],  where  each  element  in  tiie 

vector  is  given  by;  c. 

a . 

I 


i ! 

Yj  = — g— , I = 1,  2,  3,  ..., 

E'  (a.Vk!) 
k=0 


(i<.l) 


where  a.  and  c.  are  the  offered  load  and  link  capacity  of  the  ith  link, 
respectively,  t is  the  number  of  links  in  the  network. 

Clearly,  the  link  grades  of  service  are  not  invariant.  The  value 
y.  is  a nonlinear  function  of  the  offered  ioads,  a.,  and  the  number  of 
trunks,  c.,  in  each  of  the  links  in  the  netw'rk.  In  vector  notation, 

(4.1)  can  be  expressed  as: 

1 ° f.]  ('*•2) 

The  trunk  group  offered  load  vector,  in  turn,  depends  on  the  ^ vector 
and  the  traffic  distribution  vector  t.  That  is. 


The  traffic  distribution  vector  is  dependent  upon  the  traffic  matrix, 
routing  table,  call  control  rule,  and  link  blocking  probabilities.  Both 
and  are  non-linear  functions. 

Knowing  the  vector,  the  NNGOS  can  be  found  from  the  augmented 
route  trees*  for  each  node  pair.  For  example,  consider  the  augmented 
route  tree  shown  in  Fig.  4.6. 

*These  are  defined  by  the  routing  table. 
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Figure  A. 6 An  Augmented  Route  Tree 


The  NNGOS  is  found  as  follows: 

Pr{P^  used}  = 

Pr{P2  used}  = X2X^(l-x^) 

Pr{P^  used}  = Xj^x^(l-x^ ) (l-x^x^) 

Then, 

NNGOS,  „ = 1 - E Pr{P.  used} 

A-*-B  . j 

J 

= 1 - [x^  + X2X^(l-Xj)  + x^^x^(l-x^)  (l-x^x^)] 

Equation  (A. 4)  is  again  a nonlinear  algebraic  equation.  From  the 
nonlinear  relationships  between  and  node-to-node  grade  of 

services,  it  should  be  apparent  that  a two-step  algor i thm--one  step  to 
arbitrarily  assign  trunks  and  the  second  to  generate  a routing  table — 
does  not  provide  an  answer  to  the  problem  posed  in  section  3.  The  need 
for  simul taneously  generating  a routing  table  and  assigning  trunks  is  evident. 
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We  are  now  in  a position  to  answer  the  question  posed  at  the  beginning 
of  this  sub-section:  from  the  coupling  between  equations  (A. 2)  and  (A.3)f 

it  is  evident  that  information  regarding  the  desired  network  performance, 
call  control  rule  and  the  traffic  matrix  are  not  sufficient  for  the 
determination  of  the  branch  capacity  vector  and  the  routing  table  inde- 
pendently. To  determine  the  branch  capacity  vector,  the  vector  of  offered 
loads  to  the  links  must  be  known.  The  need  to  make  some  approximations 
or  assumptions  is  evident. 

The  offered  load  to  various  links  is  controlled  by  the  routing 
table  and  the  Erlang  loss  formula.  Prior  to  making  any  assumptions,  we 
note  some  of  the  properties  of  the  Erlang  B formula  [6,  10]: 

The  elements  of  the  carried  load  vector,  , can  be  obtained  as 
fol lows: 

a' . = a(l-y.)  (^.5) 

The  utilization  factor  p.  of  the  ith  trunk  group  in  steady  state 
is  defined  as: 


a 


I 


P. 

I 


('t.e) 


I 

Investigations  [6]  into  equations  (^.1),  (^.5)»  and  (I4.6)  have  led  to 
the  conclusion  that  larger  trunk  groups  are  more  efficient  than  smaller 
ones;  that  is,  as  the  number  of  trunks  in  a link  is  increased  and  the 
offered  load  is  increased  such  that  the  probability  of  blocking  remains 
constant,  the  utilization  factor  increases.  The  disadvantage  of  a link 
having  a larger  utilization  factor  is  its  higher  sensitivity  to  overload 
We  shall  use  some  of  the  above  ideas  in  devising  a scheme  for  distributi 
the  loads  to  various  trunk  groups.  Rather  than  making  some  trunk  groups 
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extremely  large  and  some  extremely  small,  we  shall  attempt  to  distribute 
the  load  such  that  the  variation  in  the  offered  loads  and  link  capacities 
is  minimal.  THp  utilization  factor  of  all  the  links  will  then  be  nearly  equal. 
By  considering  the  flow-problem,  rather  than  the  performance  analysis 
problem,  Paz  [11]  has  shown  that  improved  traffic  handling  capability  is 
obtained  if  the  variation  in  normalized  branch  flows  (ratio  of  traffic 
flow  and  capacity  of  a link)  is  minimal. 


j 4.3.2  Description  of  the  Algorithm 

? The  algoirthm  described  in  this  section  is  basically  a one-pass 

scheme  for  obtaining  the  branch  capacity  vector  and  routing  table,  given 
a BPMAX  value  and  the  traffic  demand  matrix.  Even  though  a nonlinear  analysis 
method  is  used  for  allocating  the  fully  variable  trunks,  we  shall  nevei — 
theless  make  use  of  some  linear  approximations.  Provisions  have  been 
made  to  detect  and  correct  the  d iscrepenc ies  between  the  actual  and  the 
desired  network  performance  that  may  arise  due  to  these  approximations. 

In  obtaining  this  algorithm,  the  following  assumption  [see  section 
2.4  for  others]  is  made: 

• The  time  required  to  search,  size,  and  later  release  a 
trunk  i s negl ig ibl e. 

The  algorithm  can  be  programmed  on  a digital  computer.  An 
algorithmic  flow-chart  is  shown  in  Fig.  4.7.  As  indicated  in  Fig.  4.7, 
there  are  basically  six  steps  involved. 

To  facilitate  the  understanding  of  the  various  steps  Involved,  we 
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will  use  the  seven-node  hybrid  network  of  Fig.  4.8  to  illustrate  some 
of  the  intermediate  results.  Again,  the  terrestrial  links  are  represented 
by  solid  lines.  The  satellite  ground  stations  (GS)  are  represented  by 


Assign  trunks  and  obtain  an  estimate 
of  the  network  performance. 


STEP  F 


Are  ail  NNGOS 
below  BPMAX? 


NC 

Modify  the  routing 
capacity  vector  to 

table  and  the  link 
improve  the  NNGOS. 

Figure  ^,7,  Algorithm  for  generating  a Routing  Table  and  a Link  Capacity 
vector  for  a given  BPMAX  and  Traffic  demands. 
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Figure  k.B  A Hybrid  Communication  Network 

triangles  and  the  switching  centers  (TC)  are  denoted  by  circles.  The 
fully  variable  trunks  have  not  been  shown.  Only  the  satellite  capacity 
(the  total  number  of  fully  variable  trunks  available)  is  known.  The 
traffic  matrix  and  pre-assigned  trunk  informations  are  given  in 
Table  and  (b),  respectively. 

We  will  now  make  some  definitions  and  assumptions: 

• Def ; Satellite  Section:  The  part  of  the  hybrid  system  comprising  of 

the  ground  stations  is  defined  as  the 


sate  1 1 i te  sect  ion. 
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(a)  Traffic  Matrix  In  Erlangs  for  the  Hybrid  Network  of  Fig. 


Link  Number 

Terminal  Nodes 

Number  of  Trunks 

1 

1-3 

35 

2 

2-4 

28 

3 

2-5 

25 

k 

3-4 

15 

5 

3-7 

8 

6 

3-5 

15 

7 

4-5 

13 

8 

5-7 

14 

3 

6-7 

10 

10 

5-6 

6 

Number  of  fully  variable  trunks  = 80 


(b)  Pre-assigned  Trunk  Group  Information  for  the  Hybrid  Network  of  Fig.  A. 8 


Def ; Terrestrial  Section;  The  part  of  the  hybrid  system  comprising  of 


the  switching  centers  is  defined  as  the 
terrestrial  section. 


A switching  center  and  a ground  station  in  close  proximity, 
represented  by  A,  is  considered  part  of  the  satellite  section.  Thus,  in 
Fig.  A. 8,  nodes  I and  2 comprise  the  satellite  section  and  nodes  3 through  7 
comprise  the  satellite  section. 


• Def;  Augmented  Traffic  Matrix  [t.jl  defines  the  total  traffic  between 
node  pairs  in  the  staellite  section  after  the  decision  to 
spill  traffic  to  different  nodes  has  been  made. 


* ~ originating  at  node  2 and  destined  for  node  m defined 

by  the  input  traffic  matrix  information. 


The  algorithm  for  "demand  assignment"  in  hybrid  communication 
networks  comprises  of  the  following  steps: 

STEP  A:  This  step  involves  the  computation  of  the  total  traffic  offered 

to  the  satellite  section  and  some  of  the  links  with  fixed  capacities  by 

the  terrestrial  section.  This  computation  is  easily  made  from  the 
2 

information  provided  by  the  desired  network  performance,  traffic  matrix, 
node-link  incidence  relationship  of  the  fixed  section  of  the  network. 

For  example,  consider  a switching  center  with  node  degree  exactly  one. 
Then  all  of  the  traffic  originating  at  this  node  must  be  routed  on 
this  link.  Similarly,  the  traffic  destined  to  this  node  must  also  travel 
2 

We  will  assume  the  call  control  rule  is  known.  It  is  OOC  in 
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the  same  link.  The  general  structure  of  the  augmented  route  trees  for 
such  cases  are  shown  in  Fig.  l*.?(a)  and  (b). 


|i 

; i 

ii 


I 

i 


1 


!i 


Figure  4.9  Augmented  Route  Trees  showing  (a)  traffic  away 

from  node  A (with  node  degree  1),  (b)  traffic  to 
node  A. 

Jt 

The  approximate  offered  load  to  link  connecting  A-B  in  Fig.  4.9  is 


given  by: 


"a-b  = ? ^ Vi 

I 


+ E tr. 


i-»-A 


, i / A 


(4.7) 


Once  a,  „ is  known,  the  blocking  probability  of  link  A-B  can  be  computed 
A“d 

by  Erlang's  formula.  The  additional  traffic  that  node  B must  route  is 
given  by: 


Ve<'Vi> 


&>i 


i / A (4.8) 

Is  the  additional  traffic  originating  at  node  B and  destined 


"Bh-I 
where  t 

from  node  i.  Similarly,  we  can  compute  the  additional  traffic  to  the 
other  ground  stations. 


This  approximation  is  valid  for  small  values  of  BI'MAX — \.(-,ich  is 
quite  often  the  case  wc  deal  with. 


(a)  Augmented  Traffic  Matrix  (in  Erlangs)  for  the  network  of 
Fig.  4.8  at  the  end  of  Step  A. 

1 


* denotes  a spill  switcti 

(b)  Partial  Kouting  lablo  for  the  network  of  Fig.  4.8 
at  the  end  of  Step  A. 
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In  the  case  of  SC's  with  node  degree  greater  than  one,  the  traffic 
offered  to  each  of  the  links  incident  to  it  is  in  direct  proportion  to 
their  capacities  (linear  approximation).  The  additional  traffic  that 
the  ground  stations  must  route  can  again  be  found  by  using  eqn.  (A. 8). 

When  the  decision  to  spill  the  traffic  to  the  different  ground 
stations  has  been  made,  the  appropriate  blocks  in  the  routing  table 
are  filled  with  these  choice  levels. 

To  illustrate  the  use  of  this  step,  refer  to  the  network  of  Fig.  ^.8. 
The  Hfiprox imate"*^  values  for  the  elements  of  the  augmented  traffic 
matrix  are  depicted  in  Table  The  routing  table  at  the  end  of 

this  step  is  given  in  Table  ^.5{b). 


STEP  B:  From  the  augmented  traffic  matrix,  t,  for  the  satellite  section, 

compute  the  average  load  between  a node  pair: 


1 : j 

c f c_  1 \ ^ ^ ^!->i  i ~ 1 » 2,  ...,  S 

i j ' J i = 1,  2,  ....  S 


(^.9) 


■!-"j 


tr.  . + t.  . , 

M-j  i-«-j  » 


where  S is  the  number  of  ground  stations.  Thus  for  the  network  of  Fig, 

4.8: 

<t>  = 5.74  Erlangs. 


STEP  C:  In  this  step  we  make  an  important  assumption:  the  rel iabi 1 t ies 

3 

of  all  the  GS-to-GS  links  are  equal  . Based  on  this  assumption  we  can 
compute  the  link  reliability  required  to  meet  the  prescribed  BPMAX  value. 


the  completion  of 


this  algorithm. 


these  may  turn  out  to  be  unequal. 


^These  values  are  approximate  because  some  of  the  traffic  is  lost  at  the 
switching  center  accessing  the  ground  station.  Consequently,  some  of  these 
element  values  are  larger  than  what  the  analysis  will  depict. 


^0 


We  compute  the  link  rellabllties  for  different  routing  strategies  as 
fol  lows**: 

Direct  Routing:  ^0  ^ ^ ” BPMAX 

2 

One  Alternate  Route:  + Xj(l-Xj)  = 1 - BPMAX  (i4,I0) 

2 2 

Two  Alternate  Routes:  ^ ~ ^ ” BPMAX 

Equations  (^.10)  are  nonlinear  algebriac  equations,  and  can  be 
easily  solved  by  any  fixed-point  iteration  method  [12], 

From  the  required  values  of  Xq,  x^ , and  X2*  we  can  compute  the 
average  number  of  trunks  per  link  required: 


where 


a.[  ] {1  - y.(s.,  a.)} 

I y.  (s  , a.  1 I I ’ k 

‘ k,  i = 0,  1,  2 


(A.n) 


S,  - number  of  trunks  required  in  a link  when  the  offered 
load  is  aj^  and  the  required  blocking  probability  is 
y.  (=l-x. ) . 


and 


y(s.-l,aj^)  can  be  computed  from  equation  (^.1) 


2<t> 


aj  = 2[1  + Xj  (1-x^ )]<t> 

a2  = 2[1  + X2(l-X2  )(2-X2)]<t> 


Equation  (l*.ll)  must  again  be  solved  using  the  fixed-point  iteration 
approach.  It  should  be  noted  that  s.  is  an  integer;  consequently  a proper 
error  criterion  must  be  set  in  the  Iteration  algorithm.  The  behavior  of  the 

jj 

Equations  (^.lO)  assume  that  the  alternate  paths  are  comprised  of 
two  links.  Since  we  are  dealing  with  a complete  graph  (satellite  section 
only),  this  assumption  is  quite  valid. 


Sj  vs.  NNGOS  curve  will  be  similar  to  Fig.  4.5  for  a given  traffic 
matrix.  The  choice  of  a routing  strategy  can  be  easily  made  once  we 
have  all  the  and  Information  on  the  total  number  of  variable  trunks. 

The  purpose  of  Step  C Is  to  give  an  est Imate  of  how  many  fully 
variable  trunks.  In  addition  to  the  preassigned  trunks,  will  be  required 
to  satisfy  the  BPMAX  criterion.  This  step  views  the  satellite  section  of 
the  network  to  be  symmetrical  In  routing,  traffic  matrix,  and  network 
configuration.  Such  an  estimate  Is  generally  optimistic  and,  consequently. 
In  the  actual  network,  the  number  of  additional  trunks  required  are  larger. 
But,  nevertheless,  this  step  serves  Its  purpose  by  helping  to  decide  on 
a routing  strategy. 

STEP  Pi  Once  the  decision  on  how  many  alternate  routes  will  be  used  has 
been  made  in  Step  C,  the  next  step  Is  to  form  the  routing  table  and  keep 
an  account  of  how  much  load  will  be  offered  to  the  links  connecting  two 
ground  stations. 

To  help  set  a criterion  for  forming  the  routing  table,  we  will 
make  use  of  the  results  given  In  section  4.2.  Also,  the  fact  that 
larger  capacity  links  (with  larger  offered  load)  are  more  efficient 
than  the  smaller  capacity  links  will  be  used. 


A two-step  procedure  Is  folowed  for  assigning  choice  levels  In  the 

■I 

routing  table.  These  steps  are: 

Step  1 : In  this  step,  the  first  choice  level  for  routing  of  calls 

from  ground  station-to-ground  station  is  selected.  Since  one  of  the 
configurations  available  Is  the  complete  graph,  an  obvious  candidate  for 
first  choice  level  In  the  routing  table  is  the  direct  choice.  So  the 
first  choice  for  routing  calls  between  two  ground  stations  Is  the  direct 
link.  We  fill  up  the  routing  table  accordingly. 


hi 


The  offered  load  to  the  links  connecting  every  pair  of  ground  stations  in 
the  network  is  also  recorded.  The  offered  load  to  a link  connecting  two 
ground  stations  i and  j is  given  by: 


a^|.  = t + t . 
i^J  i-M 


(^.12) 


where  the  superscript  on  a is  used  to  denote  that  it  is  the  offered 
load  due  to  path  1.  If  a direct  routing  strategy  was  chosen  in  step  C, 
step  E is  performed.  Otherwise,  we  perform  the  following  step: 

Step  2:  In  this  step  we  choose  the  second  and  third  choices  of  the 

routing  table  for  the  GS-to-GS  calls.  To  facilitate  the  assignment  of 
these  choice  levels,  we  have  established  the  following  criterion: 

The  proposed  alternate  destination  for  any  block  in  the  routing 
table  should  be  such  that  it  tends  to  minimize  the  variation  in  the 
offered  load  to  links  connecting  every  pair  of  ground  stations. 

This  criterion  can  be  easily  implemented  in  an  algorithmic  manner. 
The  offered  loads  found  in  step  1 above  can  be  arranged  in  an  ascending 
or  descending  order.  The  second  choice  levels  are  selected  such  that 
the  links  with  smaller  offered  loads  from  step  1 are  used  on  these  paths. 
The  load  offered  to  a trunk  group  due  to  the  first  alternate  path  can  be 
easily  computed.  For  example,  if  the  calls  from  node  i to  node  k use 
the  link  on  its  second  path,  then  the  additional  offered  load,  a^^. 

is  given  by: 


('*.13) 


where  x was  found  in  step  C. 


If  the  same  link  appears  on  the  second  path  of  other  node  pairs,  the 
additional  offered  loads  are  computed  by  using  eqn.  ('*.13). 


Similarly,  if  the  calls  from  node  n to  node  k use  the  link  ^ 
on  its  third  choice  level,  then  the  additional  offered  load,  's 

given  by: 


('».1M 


Thus,  in  step  2,  we  assign  second  and  third  choice  levels  (depending 
upon  which  routing  strategy  we  have  selected)  and  keep  an  account  of  the 
offered  loads  of  each  of  the  links  present  between  two  ground  stations. 


The  combined  offered  load  for  link  is  given  by; 


..  . - a E a^^,.,  • 

I -j  r^j  1'^‘j 


-P3 


IT 


('t.lb) 


where  the  summation  is  used  to  emphasize  the  fact  that  the  same  link 
can  appear  on  the  second  and/or  third  paths  of  different  node  pairs 
more  than  once;  *-’■  is  because  of  bi-directional  link  assumption. 

We  illustrate  this  step  by  once  again  referring  to  the  network  of 
Fig,  1*.8.  The  routing  table  and  the  vector  of  estimated  offered  loads 
to  the  trunk  groups  at  the  end  of  this  step  are  shown  in  Table  4.6. 


Tabic  4.6 


^\^T0 

1 

2 

3 

4 

5 

6 

7 

1 

mm 

3 

3* 

B 

3* 

3* 

2 

4* 

4 

5 

B 

3 

1 

4 

5,7 

6,7 

B 

4 

3 

2 

3,7 

5,6 

6,7 

B 

5 

3 

2 

3,7 

4,6 

6,4 

7,4 

6 

3 

5 

3,7 

B 

7,4 

7 

3 

4 

3,4 

'»,3 

5,4 

6,3 

i.e)  Routing  Tabic  at  the  end  of  Step  0 


Total  Offered  Load  = (Offered  Load  on  Path  1)  + 


X- ( 1 -x)  • (off ered  load  on  Path  ?.)  , x = .897 


^45 


STEP  E:  The  Erlang  loss  formula  was  expressed  in  vector  notation  in 
eqn.  (A, 2).  To  compute  the  link  blocking  probabilities,  the  information 
regarding  the  offered  load  and  the  number  of  trunks  in  the  link  must  be 
known . 

i*  , 

In  our  case,  the  information  regarding  the  desired  link  blocking  i 

i 

probabilities  and  the  load  offered  to  the  various  links  was  obtained  in 
steps  C and  D,  respectively.  With  this  information,  the  branch  capacity 
vector  can  be  obtained  using  equation  (it. 11).  The  number  of  fully  variable 
trunks,  in  addition  to  the  preassigned  trunks,  required  can  then  be  easily 
computed. 

Thus,  at  the  completion  of  step  E we  have  both  a routing  table  and 
a branch  capacity  vector.  There  is,  however,  a problem  associated  with 
the  design  algorithm  that  we  have  enunciated:  The  blocking  probability  of 

ground  station-to-ground  station  traffic  will  be  less  than  the  prescribed 
BPMAX  value,  but  there  is  a possibility  that  the  NNGOS  associated  with 
the  nodes  in  the  terrestrial  section  may  exceed  the  desired  upper  bound. 

For  example,  consider  the  augmented  route  tree  of  Fig.  ^.10.  The  calls 
from  node  A are  spilled  to  the  ground  station  B.  According  to  our  algorithm, 
the  trunks  are  allocated  such  that  the  NNGOS  of  the  calls  from  B to  C 
is  less  than  the  BPMAX  value.  The  blockJng  probability  of  the  calls  from 
A to  C is: 

NNGOS^  ^ = 1 - x(l-  NNG0S_.) 

A-^C  n L 

According  to  the  BPMAX  requirement,  MMG0S^_^^  ^ BPMAX*  Then, 

1 - x(l-  NNGOSg^^) 

1 - x(l-NNG0S„  ^ * BPMAX 

r 

These  are  the  links  between  a pair  of  ground  stations. 
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But  BP  * BPMAX,  and  therefore 
X “ 1 

if  the  BPMAX  criterion  is  to  be  satisfied  by  the  NNGOS  . In  general, 
X wiil  be  less  than  1 because  of  the  finite  blocking  probability  of  the 
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Figure  ^,10  An  Augment  Route  Tree 


links  between  the  switching  center  and  the  ground  station.  The  NNGOS 
of  a node  pair  using  several  tandem  switches  to  route  its  calls  may  also 
exceed  the  desired  BPMAX.  Therefore,  at  the  end  of  step  E,  the  need  to 
modify  the  branch  capacity  vector  and  routing  table  is  evident.  The 
modifications  are  carried  out  in  step  F. 

However,  before  describing  step  F,  we  will  illustrate  step  E by 
considering  the  network  of  Fig.  I|.8.  At  the  end  of  step  D,  the  routing 
table  and  estimated  offered  ioads  to  links  ^ through  13  were  obtained, 
and  are  shown  in  Table  A, 6,  From  the  offered  loads  and  the  desired  link 
blocking  probabilities  (=0.103),  we  can  obtain  the  required  number  of  trunks 
in  each  link  using  eqn.  (^.11).  These  results  are  shown  in  Table  *4.7. 
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C Table  ^.7  Trunk 

f-  f 

47 

Group  Information 

after  step  E, 

Link  Number 

Terminal  Nodes 

Numbe*"  of  Trunks 

! 1 

1-3 

Fixed  (35) 

2 

2-4 

Fixed  (28) 

J ' 

2-5 

Fixed  (25) 

4 

3-4 

25 

3-7 

13 

> 6 

3-5 

17 

1 7 

4-5 

16 

! ® 

5-7 

14 

! q 

6-7 

15 

' 10 

5-6 

17 

1 1 

4-7 

13 

12 

4-6 

10 

13 

3-6 

16 

The  analysis  based  upon  the  routing  table  given  in  Table  A. 6 and 
the  trunk  group  information  in  Table  ^,7  is  given  in  Fig.  4.11,  We  note 
the  following  salient  features  from  the  analysis: 

(1)  The  link  biocking  probabilities  of  links  4 through  13  are  not 
0.103  but  close  to  it.  This  is  to  be  expected,  since  the 
capacity  of  a link  must  be  an  integer  value. 

(2)  The  NNGOS  of  the  GS-to~GS  nodes  are  less  than  or  close  to  the 
prescribed  BPMAX  vaiue. 

(3)  Some  of  the  NNGOS  do  not  satisfy  the  BPMAX  requirement. 

To  improve  the  highest  NNGOS  in  the  network,  we  modify  the  branch 
capacity  vector  and  the  routing  table  in  the  next  step. 

STEP  F:  Before  making  any  modifications  to  the  routing  table  and/or  the 

branch  capacity  vector  obtained  in  steps  D and  E,  it  must  be  ascertained 
if  such  modifications  would  Indeed  improve  the  NNGOS  that  exceed  the 
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prescribed  BPMAX.  For  example,  the  value  of  x in  Fig.  ^.10  could  be 
much  less  than  1 so  that,  even  if  the  blocking  probability  of  calls  from 
B to  C is  made  zero,  the  blocking  probability  of  calls  from  A to  C is 
still  quite  high  ( - 1 -x) . In  such  cases , it  is  futile  to  make  any 
mod  if ication. 

But,  once  the  decision  to  make  the  modifications  is  made,  a criterion 
must  be  established  to  make  such  changes.  To  arrive  at  such  a criterion, 
we  note  that  the  number  of  trunks  in  a link  must  have  an  integer  value. 

It  is  therefore  possible  that  the  reliabilities  of  the  links  between  the 
ground  stations  may  be  higher  than  what  was  actually  required.  We  can 
insert  or  delete  choice  levels  from  the  routing  table  to  improve  the 
highest  NNGOS  in  the  network  [3].  It  is  possible  that  a few  such 
modifications  may  lower  all  the  NNGOS  below  the  specified  BPMAX  value. 

In  the  event  that  the  above  strategy  fails,  the  next  choice  is  to 
increase  the  capacities  of  the  links  appearing  on  the  paths  of  the  node 
pairs  with  the  highest  NNGOS.  Obviously,  these  links  must  be  present  between 
two  ground  stations,  so  that  fully  variable  trunks  can  be  allocated  to 
increase  their  traffic  handling  capacity. 

To  illustrate  this  step,  the  branch  capacity  vector  and  the  routing 
table  that  satisfies  the  BPMAX  requirement  (we  recall  that  the  BPMAX 
requirement  was  not  satisfied  at  tlie  end  of  step  0 is  given  in  Table  A. 8. 
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1 


Tcble  l<.8 


I 

2 

3 

it 

5 

6 

7 

1 

3* 

3 

B 

B 

3* 

2 

'♦,5 

mm 

it 

5 

5-A 

mm 

3 

1 

^5 

^,7 

5,7 

Q| 

B 

It 

3,6 

2 

3,7 

5,6 

mi 

7,6 

5 

3,6 

2 

3,7 

it,fc 

B 

B 

6 

3,7 

n 

3,7 

B 

B 

7,i» 

7 

3,6 

^,5 

3,^ 

^,3 

B 

6,3 

(a)  FIfial  Routing  Table 


Link  Number 

Terminal  Nodes 

Number  of  Trunks 

1 

1-3 

35 

t 

2-it 

28 

3 

2-5 

25 

it 

3-it 

25 

5 

3-7 

13 

6* 

3-5 

18 

7 

it-5 

16 

8a 

5-7 

15 

9 

6-7 

15 

10 

5-6 

17 

11* 

it-7 

lit 

12 

it-6 

10 

13* 

3-6 

17 

(b) 

Final  Trunk  Capacity 

Vector 

The  links  whose  capacities  were  increased  are  indicated  by  an  asterisk. 
Also,  alterations  in  the  routing  table  were  made  to  improve  the  NNf.OS, 
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An  analysis  of  the  network  based  on  the  information  in  Table  4.8  is 
shown  in  Fig.  4.12.  We  note  that  all  the  NNGOS  are  below  0.02;  thus, 
we  have  completed  our  objective  of  designing  the  network  of  Fig.  4.8 
such  that  the  BPMAX  criterion  is  satisfied. 

To  conclude  this  section,  we  will  summarize  the  algorithm  for  trunk 
assignment  and  routing  table  generation: 

The  algorithm  begins  by  assigning  choice  levels  in  the  blocks  of  routing 
table  associated  with  the  switching  centers  accessing  the  satellite  facility. 

The  traffic  spilled  to  the  ground  stations  is  recorded,  and  an  augmented 
traffic  matrix,  which  depicts  the  new  traffic  between  ground  stations,  is 
formed. 

Based  on  the  prescribed  BPMAX  value  and  the  average  traffic  (obtained  from 
the  augmented  traffic  matrix)  between  ground  stations,  the  decision  whether  to 
use  direct  or  alternate  routing  is  made.  This  choice  is  based  on  which  routing 
strategy  uses  fewer  number  of  fully  variable  trunks.  The  blocking  probabilities 
of  the  links  between  ground  stations  required  to  satisf^  the  BPMAX  specifi- 
cation (for  the  particular  routing  strategy  chosen)  are  noted. 

The  routing  table  is  now  constructed,  if  alternate  routing  is  used,  the 
alternate  routes  are  generated  such  that  the  variation  in  the  offered  loads 
to  the  links  between  ground  stations  is  minimal.  The  offered  loads  to  these 
links  are  also  recorded. 

Based  on  the  required  link  blocking  probabilities,  and  their  offered 
loads,  the  number  of  trunks  needed  in  each  link  is  computed  using  the 
Erlang's  loss  formula.  An  estimate  of  the  network  performance  is  made  to 
check  if  the  BPMAX  criterion  is  satisfied.  Further  modifications  to  the 
link  capacity  vector  and  the  routing  table  are  made  if  some  of  the  node-to- 
node  grades  of  service  exceed  the  prescribed  BPMAX  value. 
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5.  FURTHER  RESEARCH 

In  section  we  presented  a heuristic  algorithm  for  generating  a 

branch  capacity  vector  and  a routing  table  such  that  all  the  node-to-node 

grades  of  service  in  the  hybrid  network  are  below  a prescribed  value. 

Furthermore,  the  number  of  fully  variable  trunks  allocated  to  achieve 

the  network  performance  criterion  is  minimal.  To  assess  the  usefulness 

of  the  algorithm  a digital  computer  program  should  be  developed.  Then 

studies  on  some  of  the  hy  Hd  networks  in  planning  stages  can  be  made. 

Such  studies  will  show  if  this  algorithm  is  adaptable  in  near  real-time 
1 

computing  . 

The  algorithm,  at  present,  makes  a tacit  assumption  that  the  traffic 
demands  are  completely  different  everytime  this  algorithm  is  used.  In 
the  actual  operating  system,  this  assumption  may  not  be  necessary;  that 
is,  the  rate  of  call  arrivals  may  change  at  only  a few  switching  centers. 
In  such  cases,  it  is  unnecessary  to  start  the  algorithm  from  step  A; 
only  minor  alterations  in  the  branch  capacity  vector  and  routing  table 
may  suffice.  The  criteria  for  making  such  alterations  needs  to  be 
investigated.  With  a slight  modification,  this  algorithm  could  be  made 
to  detect  what  changes  in  the  traffic  matrix  have  occured,  and  what 
action  needs  to  be  taken.  Such  modifications  will  undoubtedly  reduce 
the  computation  time  Involved  in  certain  cases. 

The  assumption  that  the  network  is  in  statistical  equilibrium 
renders  the  single-moment  method  inapplicalbe  in  real-time  computation, 
where  new  traffic,  demands  are  made  every  minut>  or  so.  Consequently, 

'sased  on  the  experience  gained  from  developing  STARTUP  [3], 
it  is  envisaged  that  this  algorithm  will  be  useful  in  near  real-time 
comput ing. 
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some  other  methods  must  be  sought,  which  can  then  be  applied  to  the 
routing  problem  in  hybrid  communication  networks  for  optimum  assignment 
on  almost  every  demand  for  service. 

Previous  investigation  [3]  into  the  routing  problem  in  circuit  switched 
networks  with  originating  office  control  allowed  for  at  most  one  spill 
switch  between  a source-to-destination  pair.  Since  STARTUP  [3]  was  used 
as  an  analysis  tool  for  the  research  reported  in  this  report,  the  same 
assumption  was  tacitly  made.  With  the  availability  of  more  general  reliability 
analysis  algorithms  [13,  1^1,  study  into  the  feasibility  of  two  or  more 
spill  switches  between  node  pairs  should  be  made.  Such  an  investigation 
would  reveal  the  advantages  and  disadvantages  of  "multiple"  spill  forward 
action  in  circuit-switched  networks.  Multiple  spill  forward  action,  if 
found  advantageous,  could  then  be  used  in  the  terrestrial  section  of 
the  hybrid  network. 
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METRIC  SYSTEM 


BASE  UNITS: 

Quantity  Unit  SI  Symbol  Formula 


length 

metre 

m 

mass 

kilogram 

lime 

second 

s 

electru.  current 

ampere 

A 

thermodynamic  temperature 

kelvin 

k 

amount  of  substance 

mole 

mol 

luminous  intensity 

candela 

cd 

SUPPLEMENTARY  UNITS: 

plane  angle 

radian 

rad 

solid  angle 

steradian 

sr 

DERIVED  IINITS: 

Acceleration 

metre  per  set;ond  squared 

mis 

activity  (of  a radioactive  source) 

disintegration  per  second 

(disintegration)/s 

angular  acceleration 

radian  per  second  squared 

rad/s 

angular  velocity 

radian  per  second 

rad/s 

area 

square  metre 

m 

density 

kilogram  per  cubic  metre 

kgm 

electric  capacitance 

farad 

F 

A-s/\' 

electrical  conductance 

siemens 

S 

A/V 

electru.  field  strength 

volt  per  metn* 

V'm 

electric  inductance 

henry 

II 

V-sA 

electric  potential  difference 

volt 

V 

W'A 

elei:tric  resistance 

ohm 

V'A 

electromotive  force 

volt 

V 

W'A 

energy 

joule 

1 

N-m 

entropy 

iouie  per  kelvin 

I'K 

force 

newton 

N 

kg-mJs 

frequency 

hertz 

Uz 

(cycle)/s 

illuminance 

lux 

u 

Im/m 

luminance 

candela  per  square  metre 

cd/m 

luminous  flux 

lumen 

Im 

cd-sr 

magnetic  field  strength 

ampere  per  metre 

A/m 

magnetic  flux 

weber 

VVb 

Vs 

magnetic  flux  density 

tesla 

T 

Wb/m 

magnetomotive  force 

ampere 

A 

power 

watt 

VV 

)/s 

pressure 

pascal 

I’a 

N/m 

quantity  of  ele(  tricity 

coulomb 

C 

A-s 

quantity  of  heat 

joule 

1 

N-m 

radiant  intensity 

watt  per  steradian 

Wist 

specific  heat 

)oule  per  kilogram-kelvin 

)/kg.K 

stress 

pascal 

I’a 

N'm 

thermal  conduc  tivity 

watt  per  metre-kelvm 

W'mK 

veloi.ity 

metre  per  second 

m-s 

visi.osity.  dynami( 

pascal-second 

Pa*s 

viscosity,  kinematic 

square  metre  per  second 

m's 

voltage 

volt 

\ 

W/A 

volume 

cubic  metre 

m 

wavenumber 

reciprocal  metre 

(w8ve)/m 

work 

|oule 

1 

N*m 

SI  PREFDCES: 

Multiplu.dtion  Fac  tors 

I’rnfix 

SI  Symh 

1 000  000  000  000 

= 10’^ 

lera 

1 

1 000  mio  000 

= 10'* 

K>Hfl 

(; 

1 OOO  000 

= in'- 

mega 

M 

1 ()(M) 

10» 

kilo 

k 

100 

10* 

herlo* 

h 

10 

10’ 

deka* 

da 

0 1 

= 10  ’ 

dec:!* 

d 

0 01 

10  * 

I entt* 

( 

0 001 

I0-* 

milll 

m 

0 000  001 

10  * 

micm 

M 

0 non  000  ooi 

= 10  '' 

nano 

n 

0.000  000  000  (K)l 

10'  '* 

picn 

p 

0 (MK)  000  000  000  001 

10 

femlo 

f 

0 0(HI  fMK)  000  (HMt  000  001 

10  '• 

fltIO 

a 

To  avoided  where  poaaihli’ 


L 


